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Abstract 

We consider two variational models for transport networks, an urban 
planning and a branched transport model, in both of which there is a 
preference for networks that collect and transport lots of mass together 
rather than transporting all mass particles independently. The strength 
of this preference determines the ramification patterns and the degree of 
complexity of optimal networks. Traditionally, the models are formulated 
in very different ways, via cost functionals of the network in case of urban 
planning or via cost functionals of irrigation patterns or of mass fluxes 
in case of branched transport. We show here that actually both models 
can be described by all three types of formulations; in particular, the 
urban planning can be cast into a Eulerian (flux-based) or a Lagrangian 
(pattern-based) framework. 

Keywords: optimal transport, optimal networks, branched transport, 
irrigation, urban planning, Wasserstein distance 
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1 Introduction 

The object of this paper are two transport problems, namely the branched trans¬ 
port problem and the urban planning problem. In general terms, a transport 
problem asks how to move mass from a given initial spatial distribution to a 
specific desired final mass distribution at the lowest possible cost. Different cost 
functionals now lead to quite different optimisation problems. 

Monge’s problem. The prototype of all these problem is Monge’s problem. 
In a rather general setting, let /x+,/x_ be finite positive Borel measures on M” 
with the same total mass. The transportation of /r_|_ onto p,_ is modelled by a 
measurable map t : M" —)■ M" such that li-(B) = p,+ (t~^(B)) for all Borel sets 
B. The cost to move an infinitesimal mass d/i+(a;) in the point x to the the 
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point t{x) is given by d{x,t{x)) d^+{x)^ and the total cost is then given by the 
formula 

[ d{x,t{x))dn+{x). (1.1) 

jR" 

Usually, d{x, y) = \x — y\^ for some p > 1. 

Branched transport problem. Monge’s cost functional is linear in the trans¬ 
ported mass d/r_|_(a;) and thus does not penalise spread out particle movement. 
Each particle is allowed to travel independently of the others. This feature makes 
the Monge’s problem unable to model systems which naturally show ramifica¬ 
tions (e.g., root systems, leaves, the cardiovascular or bronchial system, etc.). 
For this reason, the branched transport problem has been introduced by Mad- 
dalena. Morel, and Solimini in |MSM03| and by Xia in |Xia03| . It involves a 
functional which forces the mass to be gathered as much as possible during the 
transportation. This is achieved using a cost which is strictly subadditive in the 
moved mass so that the marginal transport cost per mass decreases the more 
mass is transported together (cf. Figure]^. We will formally introduce branched 
transport later in Section 0 

Urban planning problem. The second problem we are interested in is the 
urban planning problem, introduced in |BB05| . Here, the measures have 

the interpretation of the population and workplace densities, respectively. In 
this case the cost depends on the public transportation network E, which is the 
object of optimisation. In fact, one part of the cost associated with a network E 
is the optimal value of (0, where the cost d depends on E and is chosen in such 
a way that transportation on the network is cheaper than outside the network. 
The other part models the cost for building and maintaining the network. A 
detailed, rigorous description is given in Section [l.2[ 

Patterns and graphs. Branched transport has been studied extensively and 
has several formulations. Maddalena, Morel, and, Solimini in |MSM03) and 
Bernot, Caselles, and Morel in jBCMOSj proposed a Lagrangian formulation 
based on irrigation patterns x that describe the position of each mass particle p 
at time t by x(p, t). The difference between both articles is that in the latter the 
particle trajectories cannot be reparameterised without changing the transport 
cost. The viewpoint introduced by Xia in |Xia03| is instead Eulerian, where only 
the flux of particles is described, discarding its dependence on the time variable 
t. A very interesting aspect of branched transport is its regularity theory as 
studied in several articles, among them |BCM08] and |Xia04| for the geometric 
regularity, |MS10| for the regularity of the tangents to the branched structure, 
[San07] and |BSlla| for the regularity of the landscape function, |BS14| for the 
fractal regularity of the minimisers. Equivalence of the different models and 
formulations are instead the topic of |MS09| . |MS13| . Branched transport can 
also be modelled with curves in the Wasserstein space as in |BBS06| . |BB10| . 
|BSllb| . 

Main result of the paper. The main result of this paper is a unified view¬ 
point of the branched transport problem and the urban planning problem. In¬ 
deed, we show that also the urban planning problem can be cast in the Eulerian 
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Figure 1: The cost per transport distance is a subadditive function of the 
transported mass m for the branched transport model = to“ for some 

a S (0,1)) as well as the urban planning model (c®’“(m) = min(am,m + e) for 
some £ > 0, a > 1). 


or flux-based and in the Lagrangian or pattern-based framework. This involves 
the consideration of new functionals which are still snbadditive in the moved 
mass, but not strictly so (cf. Figure[T|), introducing several technical difhculties. 
The main theorem, Theorem |3.3.3| proves the eqnivalence between the original, 
the Lagrangian, and the Enlerian formulation of urban planning. One advant¬ 
age of these equivalences is that now one can consider regnlarity questions in 
the most convenient formulation; as an example we show a single path property 
of optimal nrban planning networks in Proposition |3.4.f| For the sake of sym¬ 
metry we also introduce an additional formnlation of branched transport which 
is based on the transport network E as in the original urban planning formu¬ 
lation. Its eqnivalence to the pattern- and flux-based formulations is stated in 
Theorem l2.4.2l 

The following paragraphs introdnce branched transport and nrban planning 
more formally via cost fnnctionals of the transport network E. Section]^ then re¬ 
calls the Eulerian and Lagrangian formulation of branched transport and proves 
their eqnivalence to the formulation based on E. Sectionj^ puts forward the 
novel Enlerian and Lagrangian formnlation of urban planning and states their 
eqnivalence, the proof of which is deferred to the separate Section]^ Sectionj^ 
also states a regnlarity result for (a subset of all) minimisers, the single path 
property, based on the model eqnivalences. 


1.1 Branched transport 


Branched transport is the first network optimisation problem we consider. As 
already mentioned, it was introdnced in |MSM03l IXia03| to model natural and 
artificial structnres which show branching. Since it is simpler to state, we here 
take the viewpoint of network optimisation and therefore introduce a new for¬ 
mulation based on the network E. The original flnx- and pattern-based for¬ 
mulations will be introdnced in Section 2.4.2[ and their equivalence to the new 
formulation will be shown in Theorem l2.4.2l 

The network is modelled as a one-dimensional set E C M”, thought of as 
a collection of pipes through which a quantity with initial distribution /i+ is 
transported to achieve the distribution The transport cost scales linearly 
with transport distance, but is strictly subadditive in the transported mass, 
which models scale effects (i. e., transporting an amount of mass in bulk is 
cheaper than splitting it up and transporting the portions separately). Precisely, 
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the cost of moving the mass from to /r_ via a rectifiable pipe network S is 
given by 

7W“[E] = inf /c“(|F|)dH^ with c“(m) =m“. 

F:S^R"\{0} is 
div 

As it will be explained in more detail in Section the vector measure J- 
describes the mass flux, and F denotes its Radon-Nikodym derivative with 
respect to "H^LE, the one-dimensional Hausdorff measure on E. The divergence 
is taken in the distributional sense, and div identifies as a flux 

source and /r_ as a sink. The parameter a € (0,1) governs how strong the scale 
effect is, i. e., how much cost can be saved by transporting the mass in bulk. An 
optimal transport network is a minimizer of the functional A^“. 

Existence of minimisers is shown in the Lagrangian or Eulerian framework 
in |MSM03| IXiaOSl IBCM05] and many other works; regularity properties of 
minimisers are instead considered in |Xian4l IRCMnSl IRSna,j : 

typically minimisers exhibit a type of fractal regularity (see |BS14p . 

1.2 Urban planning 

The second energy functional we consider has been introduced as an urban 
planning model in |BB05| . Here, has the interpretation of a population 
density in an urban region, and /i_ represents the density of workplaces. The 
aim is to develop a public transportation network such that the population’s 
daily commute to work is most efficient. 

The transportation network is described as a collection of one-dimensional 
curves; more precisely, it is a set E C M" with finite one-dimensional Hausdorff 
measure 'H^(E). An employee can travel part of the commute via own means, 
which is associated with a cost a > 0 per distance, or via the transportation 
network E at the distance-specific cost 6 > 0 with b < a. Hence, if a travelling 
path is represented by a curve 9 : [0,1] —)■ K", its cost is given by (for ease of 
notation we here identify the path 9 with its image 0([0,1])) 

an\9\j:) + bn\9r\E). 

The minimum cost to get from a point x G K" to a point y G M" is given by the 
metric 

ds(x,y) = inf{an^(9 \i:)+ bn^(9 nS) : 9 G C^,y} (1.2) 

where 

Cx,y = {9 : [0,1] —>• K" : 9 is Lipschitz, 0(0) = x, 9(1) = y} (1.3) 

denotes Lipschitz curves connecting x and y. The Wasserstein distance induced 
by this metric describes the minimum cost to connect the population density 
/r-i_ and the workplace density /i_ and is given by 

Wdj:(M+,M-)= inf / ds(x,y)dy.(x,y). (1.4) 

It is the infimum of formula ([Tt| where we choose d(x,y) = d-s(x,y). Here, 
n(/r_|_, fj,-) denotes the set of transport plans, i. e., the set of non-negative finite 
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Borel measures on the product space K" x M" whose marginals are and /i_, 
respectively, 

n(/r+,^_) = {/r e fbm(K"' x M”) : = /r+,7r2^^ = /i_} , 

where denotes the pushforward of /i under the projection : K" x K" —>■ 

M”, {xi,X2) !-)■ Xi. 

The urban planning problem is the task to find an optimal network E with 
respect to the transport cost Wds(M+iM-) an additional penalty 'H^(E), 
the building and maintaining cost for the network. This leads to the energy 
functional 

to be minimised among all sets S C K”. Existence of minimisers has been shown 
among all closed connected E (see |BB05| or |BPSS09l Chap. 3]). Without 
requiring connectedness, existence is proved in |BPSS091 Chap. 4]. 

We will actually set 6=1 and study without loss of generality, 

since £’®’“’''(E) = 6£:f’td(i]). 

1.3 Notation and useful notions 

Let us briefly fix some frequently used basic notation. 

• Lebesgue measure. £" denotes the n-dimensional Lebesgue measure. 

• Hausdorff measure. 'H’" denotes the r-dimensional Hausdorff measure. 

• Non-negative finite Borel measures. fbm(]R") denotes the set of non¬ 
negative finite Borel measures on K". Notice that these measures are 
countably additive and also regular by |R,ud871 Thm.2.18]. The corres¬ 
ponding total variation norm is denoted by || • ||fbm- 

• (Signed or vector-valued) regular countably additive measures. 

rca(M") denotes the set of (signed or vector-valued) regular countably ad¬ 
ditive measures on K". The corresponding total variation norm is denoted 
by II • Urea- 

• Weak-* convergence. The weak-* convergence on fbm(IR") or rca(M") 
is indicated by 

• Restriction of a measure to a set. Let {X, dl, /i) be a measure space 
and Y G X with Y & A. The measure yL\_Y is the measure defined by 

AiLr(^) = /r(Any). 

• Pushforward of a measure. For a measure space (X,A4,/j), a meas¬ 
urable space (Y,Af), and a measurable map T : X ^ Y, the pushforward 
of p, under T is the measure T^p on {Y,JV) defined by 

T^p{B) = p{T-\B)) for all B GAf. 

• Continuous and smooth functions. Cc(K”) and C((°(K”) denote the set 
of continuous and smooth functions, respectively, with compact support 
on K”. 
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• Absolutely continuous functions. AC(/) denotes the set of absolutely 
continuous functions on the interval / C K. 

• Lipschitz functions. Lip(/) denotes the set of Lipschitz functions on 
the compact domain I. 

• Characteristic function of a set. Let A be a set and A C X. The 

characteristic function of the set A is defined as 

• Dirac mass. Let x G K". The Dirac mass in x is the distribution 5x 
defined by 

{6x,<p) = <p(x) for all ip G 

The Dirac distribution is the measure Sx{A) = 1 if a; S A and 5x{A) = 0 
else. 


Finally, for the reader’s convenience we compile here a list of the most im¬ 
portant symbols with references to the corresponding definitions. 

• / = [0,1]: The unit interval. 


ds, Wdj,: Urban planning transport metric and transport cost (see (1.2)- 

(pl). 


Cx^y- Lipschitz paths connecting x and y (see 

J-Q-. Flux associated with a discrete graph G (see (123). 

(F, ;B(r), Pr): Reference space of all particles (Definition |2.2.1[ ). 
X- Irrigation pattern of all particles (Definition|2.2.3[). 


[x\^Tm^{x): Solidarity class of x and its mass (Definition2.2.41. 
Irrigating and irrigated measure (Definition 2.2.6). 


r^p. Cost densities of branched transport (Definition 2.2.5) and urban 
planning (Definition |3.2.1 ). 

0: Lipschitz curves on /, 0 = Lip(/). This notation is introduced in the 
framework of transport path measures (Definition |4.1.l| . 

TPM(/r_|_,/r_): Transport path measures moving onto (Defini- 

tion|4T3 • 


2 Branched transport formulations 

In this section we will present the Eulerian or flux-based and the Lagrangian 
or pattern-based formulations of the branched transport problem and state their 
equivalence to the formulation from Section [l.l[ We begin with the Eulerian 
formulation. 
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2.1 Flux-based formulation 


We start considering the formulation given by Xia in |Xia08| . 

Let /r_|_ = J2i=i discrete finite non-negative meas¬ 

ures with Qijbj > 0, Xi,yj G M". Suppose also that they have the same mass, 

fc i 

^a, = J2bj. 

i=i j=i 

Remark 2.1.1. The object of the next definition is called transport path in 
[Xia08] . and this is the commonly used term in the branched transport lit¬ 
erature. We deliberately employ the term mass flux instead, since it does not 
only encode a path, but also the amount of mass transported. This way we 
avoid confusion when referring to actual paths as one-dimensional curves. 

Definition 2.1.2 (Discrete mass flux and cost function). A discrete mass flux 
between and p- is a weighted directed graph G with vertices V{G), straight 
edges E{G), and edge weight function w : E(G) —>■ [0, oo) such that the following 
conditions are satisfied. Denoting e~ and e'^ the initial and final point of edge 
e, we require the following mass preserving conditions, 

• Oi = J2e&E{G),e-=x, w{e) - J2e&E{G),e+=x, iOT i = 1, . . . , k, 

• bj = T,e&E{G),e+=v, “ EeeE(G),e-w{e) for j = 1, . . . , I, 

• 0 = EeeE{G),e+=v ^(<^)-EeeEiG),e-=v ^ie) for U S F (G)\{xi, ...,Xk,yi,..., yi}. 

Given a parameter a G (0,1), we define the transport cost per transport length 
as c°‘{w) = w°‘ (cf. Figure]^. The cost function Mp associated with a mass flux 
G is defined as 

M^(G) = c°‘{w{e))l{e) = Y w^(e)“ ^e), 

eeE(G) eeE(G) 

where l{e) is the length of edge e. 

In order to state the branched transport problem in the case of non-discrete 
finite Borel measures, we need to replace graphs with measures. 

Definition 2.1.3 (Graphs as vectorial measures). Gonsider a weighted oriented 
graph G. Every edge e G E{G) with direction e = can be identified with 

the vector measure Pe = ('H^Le)e, and the graph can be identified with the 
vector measure 

Tg= Y (2-1) 

eeE{G) 

All mass preserving conditions satisfied by a discrete mass flux G between 
and p_ summarise as div = p^ — p_ (in the distributional sense). 

The identification of graphs with vector measures motivates the definition 
of a sum operation between graphs that we state here for later usage. 
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Definition 2.1.4 (Sums of graphs). If Gi and G 2 are weighted oriented graphs, 
then Gi + G 2 is unique graph such that 


‘^Gi+G 2 “f ‘^G2 ■ 

Definition 2.1.5 (Continuous mass flux and cost function). Let G 

fbm(IR") of equal mass. A vector measure T G rca(K";K") is a mass flux 
between /i_|_ and /r_, if there exist sequences of discrete measures with 

M+ ^ M+) ^ M-) a sequence of vector measures Tok with ^ T, 

div —/r(l. Note that div — p,_ follows by continuity w.r.t. the 

weak-* topology. 

A sequence satisfying the previous properties is called approx¬ 

imating graph sequence, and we write (/i^, 

If is a mass flux between and pi-, the transport cost Mp is defined as 

M“(J-)= inf (hminf M“(Gfe) : (/i^,J^ gJ ^ (m+, -A)l . (2.2) 

I k—^oo J 

Problem 2.1.6 (Branched transport problem, flux formulation). Given pL+,pL- G 
fbm(]R"), the branched transport problem is 

min{Mp(J^) ; J- mass flux between pL+ and /i_}. 

Remark 2.1.7 (Existence of minimisers). A minimiser exists for /r_ G fbm(IR"') 
with compact support |Xia03| . The minimum value da(fJ,+ , P-) is a distance on 
fbm(K"), which induces the weak-* convergence (see |Xia03| L 

Remark 2.1.8. It can be shown (see |Xia03| . jBCMQl]) that a mass flux if with 
finite cost can be seen as a rectifiable set S together with a real multiplicity 
F : S —>■ (0, cxd) and an orientation e : E —)■ M", |e| = 1, such that 

F = Fe CH^LE). 

The quantity F = Fe describes the mass flux at each point in E. In that case 
we have 

M^{F) = [ F“d-H^ = [ |Frd-H^. 

Jn Jt. 

2.2 Pattern-based formulation 

In this section we recall the Lagrangian or pattern-based formulation (see |MSM03| , 
|B(1M05| . |MSI3])- 

Definition 2.2.1 (Reference space). Here we consider a separable uncountable 
metric space T endowed with the cr-algebra H(r) of its Borel sets and a positive 
finite Borel measure Pr with no atoms. We refer to (T, H(r), Pp) as the reference 
space. 

The reference space can be interpreted as the space of all particles that will 
be transported from a distribution to a distribution p_. 

Remark 2.2.2. Let {X,A4,p) and {Y,Af, v) be measure spaces. A map T : A —>■ 

Y is said to be an isomorphism of measure spaces, if 


• T is one-to-one. 
















• for every N € JV, T ^{N) G M. and /i(T ^{N)) = viN), 


• for every M G Ai, T(M) G Af and = i>{T{M)). 

Recall that if F is a complete separable metric space and Pp is a positive Borel 
measure with no atoms (hence F is uncountable), then (F,P(F),Pr) is iso¬ 
morphic to the standard space ([0,1],P([0, 1]) with m = Pr(r) 

(for a proof see |Roy88[ Prop. 12 or Thm. 16 in Sec. 5 of Chap. 15] or |Vil09l 
Chap. 1]). As a consequence, the following definitions and results are independ¬ 
ent of the particular choice of the reference space, and we may assume it to be 
the standard space without loss of generality. 

Definition 2.2.3 (Irrigation pattern). Let I = [0,1] and (F,P(F),Pr) be our 
reference space. An irrigation pattern is a measurable function y : F x / —)■ K” 
such that for almost all p G F we have Xp ^ AC(/). 

A pattern y is equivalent to x if fh® images of Pp through the maps p i—>■ 
Xp,P Xp the same. Because of that, a pattern x can be regarded as a 
map X : F —>• AC(J). 

For intuition, Xp can be viewed as the path followed by the particle p. The 
image of Xpi that is Xpil)-, is called a fibre and will frequently be identified with 
the particle p. 

Here we follow the setting recently introduced in [MS13| . 

Definition 2.2.4 (Solidarity class). For every x G K" we consider the set 

Nx = {« e r : a; G XgiJ^)} (2-3) 

of all particles flowing through x. The total mass of those particles is given by 


m^(x) =Fr([x]^). 

Definition 2.2.5 (Cost density, cost functional). For 0 < a < 1 we consider 
the following cost density, 

s^ix) = c’^{m^{x))/m^{x) = [m^ix)^-^ 


where c“ is the transport cost per transport length from Definition 2.1.2 and we 
set s^(x) = oo for m^(x) = 0. The cost functional associated with irrigation 
pattern x is 


Mp{x)= / Sa(Xp(0)IXp(0|dPr(p)dP 


(2.4) 


'rx7 


The functional Mp in the above form has been introduced by Bernot, Caselles, 
and Morel in |BCM05j . 

Definition 2.2.6 (Irrigating and irrigated measure). Let x be an irrigation 
pattern. Let iQ,if : F —>• K” be defined as Zq(p) = x{pA) and if{p) = x(p, !)■ 
The irrigating measure and the irrigated measure are defined as the pushforward 
of Pp via Zq and if, respectively, 

= M^ = (zf)^Pp. 

Problem 2.2.7 (Branched transport problem, pattern formulation). Given 
p,+ ,p,- G fbm(M"), the branched transport problem is 

min{Mp(x) : = pL+ and = ^_} . 

Remark 2.2.8 (Existence of minimisers). GivenG fbm(R”) with compact 
support. Problem 2.2.7 has a solution |MS09] . 
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2.3 Reparameterisation 


In Definition 2.2.3| one may equivalently require Xp G Lip(/) for almost all p G T, 
instead of Xp € AC(/), which is the content of Propositions 2.3.3 and 2.3.4 below. 
This becomes necessary as we will later refer to results from works using either 
one or the other formulation. In addition, it allows us to assume Lipschitz 
continuous fibres throughout the remainder of the article. 

Let us first recall the following result, whose proof can be found in [AGS081 
Lem. 1.1.4]. 


Lemma 2.3.1 (Arc-length reparameterisation for AC). Let v G AC([a, &]) and 
let L = |u(t)| dt be its length. Let 


=/a 

i{s) = infjt e [a, 6] : s{t) = s} , 


then the following holds true, 

• s G AC([a, 6]) with s{a) = 0, s{b) = L, 

• V = V o t satisfies v = v o s, v G Lip([0, L]), and lul = 1 a. e. in [0, L], 

Proposition 2.3.2 (Arc-length reparameterisation of patterns). Letx ■ Lxl -G 
K" be an irrigation pattern. Suppose x has finite cost Mp (y) < oo, and define 

s : r X / ^ [0,oo), s{p,t) = f* \x[p,t)\ dr, 

f : r X [0, oo) —>■ I U {oo}, i{p, s) = infjt G I : s{p, t) = s} , 

X : r X [0, oo) K", x{P, s) = x{P, t{P, s )), 


where for notational simplicity we define the infimum of the empty set as oo. 
Then, for almost all p G T and all s G [0, oo), y(p, •) is arc-length parameterised, 
and x(-,s) is measurable. 

The proof is similar to the one of |BCM05[ Lem. 6.2] or |BCM09l Lem. 4.1, 
Lem. 4.2]. We provide it here for completeness. 

Proof. The fact that xiPi') is arc-length parameterised for all p S T follows 
from Lemma r2.3.1l 

Since x = x° (Wr x t), its measurability properties are a consequence of the 
measurability of x and (Idr x i) and of the fact that for every null set N CT x I 
the set (Idr x is a null set in T x [0, oo). 

The measurability of Idr x t is proved as in |BCM05] and follows from the 
measurability of the map t. We now show that the set t“^(]0. A]) is measurable 
for any A S M. Let {tk}k be a dense sequence in ]0, oo). Since t is nondecreasing 
and lower semicontinuous in the variable s, we have 

OO OO 

^”^([0, A]) = Pi IJIp e T : i{p,tk) < A} x [0,4 -k b] . 

k—1 

Since (p G T : i{p,tk) < A} = (p G T : s(p. A) > tk} is measurable, we obtain 
that t“^([0, A]) is measurable, too. 
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Finally, let iV C F x / be a null set, and let i? be a Borel set such that 
N C B and (Pr ® C^){B) < 5. For almost all p € F we have 

noo pi 

/ lB(p,t(p,s))ds = / lB(p,t)|f(p,t)dt= / \B(p,t)\x(p,t)\dt. 

JQ Jo Jo 

Integrating over F, we obtain 

{Pr<S)C^){ildrxi)~'^{B))= f [ lB{p,t)\x{p,t)\dtdPr{p). 

Jr Jo 

Due to Mp{x) < oo, for every e > 0 there exists (5 > 0 such that for every set 
B with (Pr <8) C^){B) < S we have 


[ sUxiP,t))\xiP,0\^t<iPr{p) 

J B 


< e. 


Since we have that 1b{x{PiO) < 1 < ^r(r)^ “s^(x(p,t)), it follows that 


(Pr (8)/:^)((Idr X f) ^{B))= j J 1b(p, t)|x(p, t)| dtdPr(p) 

< Pr(r)i"“ [ sl{xiP,t))\x{p,t)\^tdPr{p) < Pr(r)^ 

Jb 


e. 


Choosing e arbitrarily small gives (Pr(8)/l^)((Idr x i) ^{N)) = 0 as desired. □ 


We may further reparameterise the irrigation pattern. 

Proposition 2.3.3 (Constant speed reparameterisation of patterns). Let x ■ 
Fx/ —>■ M" be an irrigation pattern with finite cost Mfi{x)i Kp) — fo \xiP^ 0\ dt 
be its fibre length, and let x be as in Proposition \ 2.3.^ Then x : T x / —>■ K", 
(p,s) I—)■ x{p,s/l{p)), is an irrigation pattern which reparameterises the fibres of 
X with Xp € Lip(/) and constant velocity |xp| for almost all p € F. 

Proof. This follows from the properties of x- □ 

Proposition 2.3.4 (Reparameterised patterns have the same cost). Let x ■ 

F X / —> K" be an irrigation pattern with finite cost Mp{x) ond let x be its 
Lipschitz reparameterisation. Then Mpfx) = Mp(x). 


Proof. The proof is straightforward, once one notices that the solidarity classes 
(2.3) do not depend on the parameterisation. □ 


2.4 Equivalence between the formulations 

It has been proved by Bernot, Caselles, and Morel in |BCM08l Sec. 6] that the 
pattern-based formulation is equivalent to the formulation by Xia, even though 
Xia’s formulation does not include the particle motion, while in the pattern- 
based formulation by Maddalena, Morel, and Solimini the speed of particles 
occurs in the functional. In particular, minimisers exist for both models, and 
they can be identified with each other. 
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Definition 2.4.1 (Branched transport energies). Given two measures S 

fbm(]R") of equal mass, for an irrigation pattern x, a mass flux T, and a recti¬ 
fiable set S C R" we define 


[x] = (x), [J-] = (.F), 


where My{IF) and Mp(x) are given by (2.2| and (2.4|, respectively, as well as 


Al“[x] if M+ = F-i- and /i^ = 

cxD else, 

if dlyF" = 
oo else, 

= inf{M“’^+’'"”[J^] : F- = F'H1 lE, F : S 




\{ 0 }}. 


The last functional corresponds to the new formulation of Section p~T| Note that, 
if S is not rectifiable, then = oo (see |Xia04l Proposition4.4]). 

Theorem 2.4.2 (Equivalence of branched transport energies). The minimisa¬ 
tion problems associated with Definition \2.4-l\ are equivalent in the sense that 

minM“’'"+’'"- [x] = minM“’'"+’'"- [F] = minM“’'"+’'"- [E] . 

X s 

The optima can be identified with each other via 

Yj = {x G R” : rn^(x) > 0} , F = FH^L-Y for the density F = m^e , 
where e is the tangent unit vector to Y. Moreover, 


ip ■ dF = 


J j fiXpit)) ■ Xp(0dtdPr(p) for all p G Cc 


Proof. The equivalence of the pattern-based formulation min^j. [x] to 

Xia’s formulation minjr [F] has been proved by Bernot, Caselles, and 

Morel l IBCMOSl Sec. 6] or [BCMOQl Chap. 9]). Furthermore, for an optimal X) 
the set 

E = {a; G R” : m^{x) > 0} C jj XpW 

pGT 

is rectifiable |BCM05l Lem. 6.3], and thus for T-L^-a. e. point x h as a ta ngent unit 
vector e(a;). Defining a multiplicity via F = (see Definition 2.2.4) we obtain 
a flux F = FeH^\_Y as in Remar k [2.1.8[ and the proof of [BCMOQ] Prop. 9.8] 
implies that this flux is optimal. 

The equality minjr ]F] = mins [E] follows by choosing 

E as the rectifiable set from Remar k [2. 1 .8] corresponding to the optimal F. 

Finally, using the relation between the optimal F and x, for p G Cc(R”; R") 
we have 


p ■ dF = 


Xp(G) 


p{x) ■ F(x) dH^(x) 


x]^p(x) ■ e{x)d'H^{x) = J^p{xpit)) ■ Xp{t) dtdPr{p). 


This formula follows noting that if two fibres Xp and Xq coincide in an interval, 
then their tangents coincide TL^-a.e., too. □ 
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2.5 Regularity properties 

Due to proof of equivalence one can examine regularity properties of minimisers 
in the most convenient formulation. The following is based on patterns. 

Definition 2.5.1 (Loop-free paths and patterns). Let d G Lip(/) and let y be 
an irrigation pattern. Following |MS13I Def. 4.5], we say that 6 has a loop if 
there exist ti <t 2 < such that 

0{ti) = eits) = X , 0{t2)y^X-, 

else we say that 0 is loop-free, x is said to be loop-free if Xp i® loop-free for 
almost all p € r. 

Definition 2.5.2 (Single path property). Let y be a loop-free irrigation pattern 
and let 

= ■ Xp^{x) <Xp^{y)} ■ 

Following |BCM08| Def. 3.3] and |BCM09l Def. 7.3], x has the single path prop¬ 
erty if for every x,y with 7r(r^) > 0, the sets x(p, [Xp^{x),Xp^{y)]) coincide 
for almost all p G F^. 

Note that under the single path property, almost all trajectories from x to y 
coincide, but they need not coincide as functions of time (since the time variable 
can be reparameterised). 

Remark 2.5.3. Optimal patterns are loop-free and enjoy the single path property 
(see |BOMn8l Sec. 3], |BOMn9l Chap. 4] or |Mm3l Thm.4.1]). 

3 Urban planning formulations 

Here we will employ the same notions as in the previous section to provide 
the Eulerian or flux-based and the Lagrangian or pattern-based formulations of 
urban transport. These will then be proved equivalent to the original definition, 
e. g. from |BPSS09] . 

3.1 Flux-based formulation 

Let G = {V(G), E(G), w) be a discrete mass flux between discrete measures 
G fbm(IR"). Let E be a subgraph of G; E is not required to be connected. 
Given parameters e > 0, a > 1, the cost function Ey°‘ is defined as 

F;p’“(G,E)= ^ aw{e)l{e)+ ^ {w{e) + £)l{e), 

e£E(G)\E(T.) e£E{Y.) 

where l{e) is the length of edge e. Ep‘^{G, E) is the cost for employees to travel 
from an initial distribution of homes to a distribution of workplaces via 
the network G using public transport on E. We wish to minimise ifp“(G, E) 
among admissible pairs (G, E). For a pair to be optimal one must have 

• aw{e) < w{e) -I- e if e G E{G) \ G(E), since otherwise the pair (G, E U {e}) 
has a lower cost, and 

• aw{e) > w{e) -I- e if e G E{E), since else (G, E \ {e}) has a lower cost. 
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As a result, the cost of an edge e G E{G) for an optimal (G, S) is given by 
min(aiy(e), w(e) +e) so that the problem can be the restated with just the mass 
flux variable G. 


Definition 3.1.1 (Cost function, flux formulation). Let G = {V{G),E{G),w) 
be a discrete mass flux between discrete measures G fbm(M"). Given 

parameters £ > 0, a > 1, we define the transport cost per transport length as 
E’°‘{w) = min(ar/;, w + £) (cf. Figure]^. The cost function Ep°' associated with 
a mass flux G is defined as 

£'p“(G) = c^’°‘{w{e))l{e) = m.m{aw{e),w{e) + e)l{e), 

eeE(G) eeE{G) 

where l{e) is the length of edge e. 

If G rca(]R") is a general mass flux between general measures G 

fbm(K"), the cost function is defined as 


=inf{liminf£;;’“(Gfc) 

fc—>-oo 




Problem 3.1.2 (Urban planning problem, flux formulation). Given G 

fbm(]R"), the urban planning problem is 


min{ifp“(J^) : J- mass flux between /r_|_ and /r_}. 


Remark 3.1.3 (Existence of minimisers). The existence of mass fluxes with fi¬ 
nite cost follows from the existence of irrigation patterns with finite cost (Re¬ 
mark 3.2.31 and Proposition |4.3.l| later. Furthermore, Ep°‘ is weakly-* lower 
semicontinuous by definition, and it is bounded below by || • ||rca (since it is the 
relaxation of a functional, defined only on discrete mass fluxes, which satisfies 
the same property). Thus, graphs with uniformly bounded energy are weakly-* 
precompact, and existence of minimisers follows via the direct method of the 
calculus of variations. 


Remark 3.1.4. Note that, just like c“ for branched transport, the function 
c'^’“ is subadditive (cf. Figure]^, since a concave function whose graph passes 
through the origin is subadditive. This leads to an economy of scales and thus 
to branched structures. However, unlike c“ it is not strictly subadditive, so 
there is a slightly weaker preference for branching structures. In particular, the 
minimisers need not be finite graphs away from the support of the initial and 
final measure, and mass fluxes can locally be absolutely continuous with respect 
to Lebesgue measure £" (see Figure]^. 

Remark 3.1.5. For finite graphs, the corresponding optimal network subgraph 
E is the graph whose edges are 


E{Yf) = {e G E{G) : aw{e) > w(e) -I- e} . 


3.2 Pattern-based formulation 

Definition 3.2.1 (Cost function, pattern formulation). Let (r,H(r),Pr) be 
the reference space and let x : F x [0,1] —> M” be an irrigation pattern. For 
e > 0 and a > 1, consider the density 

re,a(x) = c^’“(TOx(a;))/mx(a;) = 


14 


if rn^(x) > 0, 
if rn^(x) = 0. 







Figure 2: Sketch of an optimal urban planning mass flux which is absolutely 
continuous with respect to Lebesgue measure in some regions. The grey shade 
indicates the local flux density. 


The cost functional ifp “ is 

'Ep“(x) = [ rl^{Xp{t))\Xp{t)\ dPr(p) dt. 

Jrxi 

Problem 3.2.2 (Urban planning problem, pattern formulation). Given /r_ S 
fbm(K"), the urban planning problem is 

min{ii'p“(x) : = /r+ and = /r_} . 

Remark 3.2.3 (Existence of a finite cost pattern). An irrigation pattern with fi¬ 
nite cost Up “ for a given pair /r_|_, of finite Borel measnres with same mass and 
bonnded snpport can readily be constrncted based on the Monge-Kantorovich 
problem. Indeed, it is well-known that the 1-Wasserstein distance between 
and is bounded (see e.g. |Vil03l Chap. 1]), 


Wi(M+,Ai-)= inf 


\x - y\dn{x,y) < oo , 


and that the inflmnm is achieved by a minimising measnre y, € n(/r_|_,/r_), 
where n(/r+, /i_) denotes the set of transport plans as in Section 1.2 


By |Roy88[ Prop. 12 in Sec. 5 of Chap. 15] there exist a measnre v on [0,1] and 
an isomorphism (p : ([0,1], R([0,1]), ^) —>■ (M" x K",R(K" x K"),/r) of measnre 
spaces. Define ip : [0, to] —>■ ]0,1] as the pseudo-inverse of the cumnlative func¬ 
tion oi v. It is clear that the pushforward of the Lebesgue measure under (po if 
is y. Now take the reference space (r,6(r),Pr) = ([0, to], R(]0, to]), £^L]0, to]) 
and define the irrigation pattern 


xivG) = Ct{p{ip{p))) with Ct:M”x 


Ct{x,y) =ty + {l-t)x. 


Since Cq and Ci are the projection on the first and second argnment, it is a 
straightforward exercise to verify that yf^ = /i+, yfi = y-. Moreover, we have 


£^p“(x) = / r^,aix{pG))\x{P,t)\dPr{p)dt < 
Jrxi 

= [ a\Ci{p{if{p)))-Co{pi'ifip)))\dPr{p) = [ 

Jr 


a\xiP,t)\dPrip) dt 


fFxI 


a\x-y\dy{x,y) = aWi{y+,y-). 

R" 
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Remark 3.2.4 (Existence of minimisers). The existence of patterns with minimal 
urban planning cost will follow from Remark |4.2.3| via the equivalence of different 
energy functionals, one of which admits a minimiser. 

Before considering the equivalence between the different formulations, let us 
state a few properties of the cost functional for later use. 


Proposition 3.2.5 (Constant speed reparameterisation of patterns). Irrigation 
patterns of finite cost can be reparameterised such that Xp G Lip(^) |xp| is 
constant for almost all p without changing the cost ifp^. 


Proof. The proof is analogous to the proofs of Propositions |2.3.^ to |2.3.4 merely 
replacing the estimate 1 < Pr(r)^““s^(x(p,t)) by 1 < rf aixiPji))- ^ 


The below closely follows |MS091 Lemma 4.4, Lemma 4.5]. 


Definition 3.2.6 (Pointwise convergence). Let Xn be a sequence of irrigation 
patterns. We say that Xn pointwise converges to x if for almost all p G T the 
curve Xn(p, ■) converges uniformly to x- 


Proposition 3.2.7 {m^ is u.s.c. and is l.s.c.). Let Xn be a sequence of 
irrigation patterns such that Xn ^ X pointwise. Let tn & I such that —>■ t. 
Then, for almost all p € T, 


m^ixiP^t)) > linisupmx„(Xn(p,tn)), (3.1) 

n—^oo 

r^aixiP,*)) < liniinfr^-(x„(p,t«)) • (3.2) 

Proof. Fix p £ r such that Xp G AC(/) (we only discard a null set of fibres) and 
define the sets A and An as 

A = [^ , An = IJ [Xfc(Pi ffc)]xfc ■ 

n k>n 


Recall that A = limsup„ A„, Pr([Xk(p,tk)jxk) = ^xk(Xk(p,tk)) and 
Pr(A) =limPr(A„) > limsupm^^(xn(p,tn)) ■ 

^ n 


We want to show that Pr(A\ [x{pA)]x) — b) that is, A C [x{p,t)]x to a 
negligible set of fibres so that m^{x{p,t)) > Pr{A) > limsup„ (xn(p, tn)), 
proving inequality (3.11. Let then q G A such that Xn{q, ■) —>■ x)?) ’) uniformly 
(we only discard a null set of fibres of A) . Recall that q G A if and only if there 
exists an increasing sequence of integers rife and a sequence Sn^, G I such that 
Xn^TSuk) = XnkiPAnJ- Supposc now by contraposition that q ^ [x(p,t)]x- 
Because of the continuity of x we have d = dist(x(p, t),xiQ, I)) > 0. Because of 
the fact that Xn X uniformly, for large k we have also dist(x(p, t), Xn,, {q, I)) > 
|, contradicting Xn^ (q, SnJ = Xn^ (P, tnj ^ x{P, t)- _ 


Inequality (3.21 follows immediately from inequality (3.1) and the definition 


of r^ 

' €,a’ 


□ 


Proposition 3.2.8 (Lower semicontinuity of Ep°‘). The functional ifp“ is 
lower semicontinuous with respect to pointwise convergence of patterns. 
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Proof. Let Xn be a sequence of irrigation patterns converging pointwise to the 
pattern x- For a given integer n and a given fibre p, define /i„ = |Xn(_p, •)! dt 
and p = \x(p, ■)\dt. As a consequence of the uniform convergence of Xn we have 
Xnip, •) -X xiPi ■) ill tbe distributional sense and thus 

< liminf/i„(A) 

n—^oo 

for any open A C I. Thanks to |MS18[ Def. C.l and Thm. C.l] and Proposi- 
tion l3.2.7l we thus have 

Ji Ji 

Integrating with respect to and applying Fatou’s Lemma ends the proof. □ 


3.3 Equivalence between the formulations 

Definition 3.3.1 (Urban planning energies). Given p,+ , /r_ G fbm(IR") of equal 
mass, for an irrigation pattern x, a- mass flux P', and a set S C K" we define 




as well as 




I i£’®’“[x] if P+ = M+ and = P-, 

I oo else, 

\£‘^'°‘[P] \i dwP = p+— P-, 
l(X) else, 

\xXdj^[p+,P-) + if S is rectifiable, 

1 00 else. 


with Wdj, defined in Section 1.2 


Remark 3.3.2. Recall that [S] has a minimiser, thanks to |BPSS091 

Lem. 4.10, Lem. 4.11, Prop. 4.15, and Thm. 4.26]. 

The next theorem is the main result of this paper. Its proof will be the 
object of the next section (Section . 


Theorem 3.3.3 (Equivalence of urban planning energies). The minimisation 
problems in Definition \3.3.1\ are equivalent in the sense that, for p- of equal 
mass and with bounded support, they possess minimisers and satisfy 


minf^.^./^+.M- r^i = min [P] = min£:®’“’'"+’'"- [E]. 

X . 7 ^ S 

Similarly to branehed transport, there are optima x> ^> OLnd S that can be iden¬ 
tified with each other via 


[ p-dP=[ [ ip{xp{t)) ■ Xp{t)^t^Pr{p) for all ip G Cc 
JR" Jr Ji 


S = {a; € K" : rn^{x) > -rj} . 
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(3.3) 

(3.4) 











TOi + m2 



mi + m2 


frii+m 


mi 


mi 


Figure 3: Optimal particle motions for the irrigation problem from Remark |3.4.2| 
(filled circles denote sources, open ones sinks, each arrow is labelled with its mass 
flux). Any 0 < m < m 2 yields an optimal mass flux or irrigation pattern. 


3.4 Regularity properties 

A consequence of the above-stated equivalence between the different models is 
the fact that regularity issues can now be considered in the most convenient 
formulation. As an example, we state the following single path property of 
minimisers to the urban planning problem. Its proof is given at the end of 
Section|421 


Proposition 3.4.1 (Single path property for the urban planning problem). 
There exists an optimal irrigation pattern x for which has the single 

path property (see Definition 2.5.2). 


Remark 3.4.2. Unlike for branched transport there exist optimal irrigation pat¬ 
terns not satisfying the single path property, for instance for 


with mi > > m 2 and s = Va^ — 1 (see Figurej^. 


4 Proof of Theorem |3.3.3| and of Proposition |3.4.1| 


In this section Theorem |3.3.3 is proved in four steps, corresponding to Propos- 
itions|4)2)T] |4.2.2[ |4.3.1 and 4.3.6 The proof of |3.4.I| follows as a corollary of 
those of Propositions 4.2. 1| and 4.2.2[ We first introduce some necessary notions 
concerning measures on the set of paths and cycles in discrete graphs. 


4.1 Transport path measures and other preliminary defin¬ 
itions and results 

Definition 4.1.1 (Transport path measures). Let 0 = Lip(/) be the set of 
Lipschitz curves I —)■ K" with the metric 

c?e( 0 i, 02 ) = inf ■|max|0i(t) — 02(<p(t))l : tp : I ^ I increasing and bijectivej’ . 

Following |BPSS09l Def. 2.5], a transport path measure is a measure 77 on 0 (en¬ 
dowed with the Borel algebra). If G fbm(]R"), we say that the transport 

path measure 77 moves onto /i_ if 

Po^V = p+, pi^r] = p-, 

where, given t G I, pt '■ Q ^ R" is defined by pt(P) = Oft). We denote by 
TPM( 7 i_|_,/r_) the set of transport path measures moving 71 + onto /r_. 
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Remark 4.1.2 (Compact sets in 0). The following compactness result can be 
obtained via the Ascoli-Arzela Theorem (see |BPSS09l p.7]). Let 0i, 02 , • ■ • be 
a sequence in 0. Suppose that the 0„ have uniformly bounded lengths and 
0n(O) G for a compact subset Q. of K". Then, the sequence On is relatively 
compact in 0 . 

Definition 4.1.3 (Parameterisation of transport path measures). Given a ref¬ 
erence space (r,yB(r), Pr), a parameterisation of a transport path measure rj is 
a function y : T —0 such that X#Pr = ??• With little abuse of notation, we 
write x(p, t) instead of x(p)(t) (the position of particle p at time t). 

Notice that the map x : T x J —>■ K" satisfies 

• x(P, •) € 0 for a. e. p e T, 

• x(’>0 i® measurable for all t G I. 

Remark 4.1.4. Note that a parameterisation of a transport path measure al¬ 
ways exists thanks to Skorokhod’s Theorem (see |BCM09( App. A, Thm. A.3] or 
|Dud021 Thm. 11.7.2]). 

Definition 4.1.5 (Cost of a transport path measure). Let a > 1 and E be a 
Borel set with 77^(S) < oo. Following |BPSS09l Chap. 2, eq. (2.2)], we define 
the cost functional as 


C^iv)= / api(0(/)\E)-bpi(0(/)nE)d77(0). 

J0 

The following is a slight refinement of |BPSS09l Cor. 2.12 and Prop. 2.14]; 
we thus only repeat the relevant parts of the proof. 

Proposition 4.1.6 (Optimal transport path measures). Given a bounded Borel 
set E C M" with P^(E) < oo and G fbm(]R") with equal mass and 

compact support, there exists a minimiser rj G TPM(/i+,/i_) of such that 

rj-a. e. path 0 G & is loop-free. Furthermore, 

C^fq) = Wd 5 ,(/r+,p_) . 

Proof. Using the notation of |BPSS09I Prop. 2.14], let us set (5e(0) = a'H^{0{I)\ 
E) -I- n E). By |BPSS09I Prop. 2.18] it is easy to see that equals its 

relaxation so that 

cm)= [ S^{9)dm= [ fe(0)d7)(0)=:C|(7)) 

Je Je 

for all transport path measures fj. By |BPSS09l Prop. 2.14], and thus also 
eg possess a minimiser q G TPM(^+, ^_). 

Next, employing the notation from Section [l.2| for a transport plan 7 G 
n(/r+,/r-) define /^(x) = d^{x,y) d'j{x,y). By |BPSS09i Prop. 2.14], 

there exist a (minimising) transport plan 7 and a (minimising) transport path 
measure q such that 

Wd,(/r+,Ar-) =/e(7) = ^( 0 ) = Cg(7). 
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Finally, let us show that ij can be chosen such that 77 -a. e. path is loop-free. 
To this end, let il C M" be a closed ball whose interior contains the support of 
and /i+ as well as E, and define 

&(0) = + l)lR.\s(0(t))) \m\dt. 

Jo 


Note that (5 e : 0 —?> [0, 00 ) is lower semicontinuous. Thus, by the same proof as 
for |BPSS09l Cor. 2.11], for any given x,y G a minimiser of Je exists in 
Cx,y = {0 G Q : 6*(0) = x,0{l) = y} with 9x,y{I) C fl. Therefore, by exactly 
the same proof as for |BPSS09l Cor. 2.12] there is a Borel function 5 : fl x fl —>■ 0 
such that 

^ 2 ( 9 ( 2 ;, 2 /)) = nun 5^{e). 

C' G Oa; ^ y 

Now assume there are x,y G ft such that q{x, y) has a loop, that is, q{x, y){ti) = 
q(x, y){t 3 ) = z and q{x, y){t 2 ) ^ 2 : for some 0 < < ^2 < ^3 < 1- Then 


S^iqix,y)) > S^{0) for e{t) 


Q{x,y){t) if t e [0,ti] U [fa, 1], 
z else. 


which contradicts the optimality of q{x,y). Thus, q(x,y) is loop-free and there¬ 
fore injective up to reparameterisation, and we have 


h{qix,y))= [ {1 +{a-l)l^r^\^{q{x,y){t)))\q{x,y){t)\dt 
Jo 

= y)(/)) -f (a - l)n^{q{x, y){I) \ E) = 6s{q{x, y)). 

If we can show 5^{q{x^ y)) = minggc^ ^s(9), then, letting 7 G n(p,_|_, /i_) be an 
optimal transport plan, as in the proof of |BPSS09l Prop. 2.14] it follows that 
r] = q^"/ is an optimal transport path measure. 

We close by proving 6s{q{x,y)) < minggc^ (the opposite inequality 

holds trivially). By |BPSS09I Cor. 2.11], 6^ = possesses a minimiser 9 G C^^y 
If 9 has a loop, that is, 9{ti) = 9{t3) = z and ^(^ 2 ) ^ z for some 0 < ti < t 2 < 
fa < 1, and if ^(^ 2 ) ^ 9[I\ [^ 1 ,^ 3 ]), then 


5^{9)>5^{9) for 9{t) 


9{t) if f G (0,ti] U (fa, 1] 
z else. 


which contradicts the optimality of 9. Thus, 9{I) must be homeomorphic to / 
and can be parameterised by an injective 9x^y G Therefore we have 

nun 5^{9) = 5y.{ 9) = 6^{9x,y) = S^i9x,y) > S^{q{x,y)) 

t't Ox j y 


as desired. □ 

Finally, let us consider the relation between discrete graphs and transport 
path measures. 


20 










Definition 4.1.7 (Mass flux cycle). Let G be a discrete mass flux. A cycle C 
is a collection of directed edges {ei,... ,ek} C E{G) such that ei U ... U Cfc is 
homeomorphic to a circle with constant orientation. 

The weight of the cycle is defined as wc{C) = rc(ei). Let = 

{e S G : wo{e) = wg{C)} be the set of edges with minimal weight in G. 

The G -reduced mass flux Gc is the discrete mass flux such that the set of its 
edges is E{Gc) = E{G) \ Jq with weights wcd^) = waie) — wg{G) for e G G 
and WGcis) = wcie) else. Notice that divGc = divG so that initial and final 
measures are unchanged. 

The cycle-reduced mass flux is the mass flux G reduced by all cycles. 

Remark 4.1.8. The cycle-reduced mass flux is well-defined since a discrete mass 
flux has at most finitely many cycles which can be reduced one-by-one. In doing 
so, it is easy to see that the reduction order does not matter. 

Lemma 4.1.9 (Cost of mass flux cycle). Let G be a discrete mass flux between 
G fbm(M") and G the cycle-reduced mass flux. Then 

(G) < (G) - ||.Fg --^Cllrca . 

Proof. Let G be a cycle of G. We have 

\\Egc - J^cWrce. = ^ WG{e)l{e) + {wG{e) - WGc{e))l{e) 

eGjg e^E(Gc) 

^ Y '^G{C)l{e) + Y WGiG)l{e) = WGiC)YKe), 

eGJg eeE(Gc)nC eGC 

where l{e) denotes the length of edge e. Likewise, since c®’“(ui) — — wq) > 

Wo for all w > Wo > 0, 

A^’“(G) - £;^’“(Gc) 

= ^ c''’“(w;G(e))/(e) + Y (c®’“(w^G(e)) - c''’“(wGc(e)))^(e) 
eGJg e(^E(Gc) 

= ^ c''’“(w;G(G))/(e) + Y {c^’‘"{wG{e)) - c^’‘^{wG{e) - WG{G)))l{e) 
e&J§ e(iE{Gc)r\C 

>WG{G)YKe) = \\EGc-m.o.. 

eec 

The result now follows repeating this procedure over all cycles and using the 
additivity of || • ||rca with respect to cycle removal. □ 

Remark 4.1.10. By |AM093l Thm. 3.5 and Prop. 3.6] or |GDL14I Thm. 1], a 
discrete mass flux between G fbm(K"’) without cycles can be identified 

with a transport path measure 77 moving 77 + to /i_. 
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4.2 Equivalence of the pattern- and set-based formulation 

The combination of the following two propositions proves the existence of an 
optimal pattern x a-nd 


X s 


as well as relation (3.41 for at least one pair of minimisers, as detailed in Re- 
mark lTT^ 


Proposition 4.2.1. For any S C K" there exists an irrigation pattern x such 
that 

< ge,a,n+,U- p]. 

Proof. By Proposition |4. 1.6] there exists an optimal transport path measure ij € 
TPM(/x+,/r_) such that rj-a.e. path 0 G 0 is loop-free and 


Wd^(/r+,/i_) = f an\e{i)\i:) + n^{e{i)ni:)dr,{9). 

Je 

Let X be a parameterisation of the optimal transport path measure 77 so that 
V = X#Pr- 

First we derive the formula 


= 


1 


lrJ{tei : x(p,t)GS:} ™x(x(P;0) 

where we introduced E C S as 

E = {x G E : m^(x) > 0 } . 


\x{p,t)\dt dPr(p), 


Since x is loop-free (thanks to Proposition 4.1.6 1 and thus Xp is injective (up to 
reparameterisation) for a. e. p G P, we have 


'{tG/: x(pG)GS} ''^x(x(P)^)) 


lx(P.i)|dt 


= I,„ns = L ^ 

Using this as well as the identity ^ /p dPr(p) = we 

obtain 


1 

rJ{tei : x(p,t)et} ^xixiPit)) 

1 


|x(p,t)|dtdPr(p) 




/r "‘x 


(x) 


lxp(/)nf:(^)‘i^r(p) d'H^(x) ='Hi(E) 


1 _ 1 


(x) ■fr ^xp{r)nE(^)d'Pr(p) - ~ 
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after application of the Fubini-Tonelli Theorem. Next we notice that 

[ IXpWldi, 'H^(Xp(/)nl]) = / \Xp{t)\<it 

{xpit)i^} "'{xp(t)es} 

for a. e. p € r due to the injectivity of Xp so that in summary, the urban planning 
cost can be estimated as 


^ / an\9{i)\^) + n\0{i)nY.)d'n{0) + en\^) 
Je 

= ^ an\xp{I) \ S) + n\xp{I) n s) dPr(p) + eP'(S) 
> ^ Cin\xp{i) \ S) + n\xp{i) n s) dPr(p) + en\^) 




IXpWIdP 


+ 


= / a / 

dr J{xp{t)m 


/{Xp(i)6S\E} 
£ 


{xpWes} m^iXpit)) 

IXpWldt 


IXpWldf 

|Xp(t)|dt dPr(p) 


{xp(t)eE} 


IXpWldi 


'{xp(t)es} 


1 + 


^x(Xp(0) 


IXpWldt 


/{Xp(i)6S\S} 


|Xp(t)|dt-a 


/{Xp(t)eE\E} 


|Xp(t)|dtdPr(p) 


> 


fVxl 


r?,a(XpW)IXpW|dPr(p)dt + f f {1 - a)\xp{t)\dtdPr{p) 

dr Jtx„(Oes\E> 


> 


Jrxi 

Thus the pattern x satisfies 


'r J{xpWes\s} 

r?.a(Xp(t))IXp(t)l dPr(p) dt = . (4,1) 


[^] < ge,a,p+,p_ p] ^ mm 5'^’“’^+[S] , 

concluding the proof. □ 

Proposition 4.2.2. For any irrigation pattern x> the rectifiable set E = {a; € 
M" : rn^{x) > so,tisfies 

ge,a,p+,p. p] _ 

Proof. Let x be an irrigation pattern. Let us define rj = X#Pt (x here is a map 
r —> 0). By definition of E, we have 


r^,a{x) 



if a; G E, 
otherwise 


and thus 

/^?.a(XpW)IXpW|dt= / (^+ m Mpt)) ) \xip,t)\<^t+ [ a\x{p,t)\dt. 

d{xp(t)es} ^ hxvidm 
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First, notice that 


/ / \x{P,t)\dt+ / a\x{p,t)\dt\ dPr{p) 

Jr \J{xpm^} J{xp{t)m J 

> [ aHi(0(/)\S)+Hi(0(/)nS)d77(0) > Wd^(M+,M-) 

Je 

by Proposition |4.1.6| Furthermore, 

[ f dtdPr(p) > f [ —^d'H^(a;)dPr(p) 

Jr i{xp(t)es} m^iXpit)) Jr JxpiPns m^{x) 

^dPr(p)d'H^(a;) = f dH^ix) = 'H^(S) 

.(x) 7s Wx(x) 


s J{per : xexp(i)} "^xV 


so that 


£’e:,a,//+ ,/x_ 


[x] = ^^?'?.a(XpW)IXpW|dtdPr(p) 

> Wd, (a^+, A^_) + eH'(E) = [E], 


concluding the proof. 


□ 


Remark 4.2.3. Due to the constructive nature of the above proofs, the existence 
of an optimal E by [BPSSOQl Sec. 4.4 and Thm. 4.26] implies the existence of 
an optimal irrigation pattern x- Indeed, let xs be the pattern constructed in 
Proposition |4 .2.1] from an optimal E, then 

inf ’'"-[x] < £:®’“’'"+’'"-[xs] < min ’^-[E] < inf M , 

X Ex 


where the last step follows from Proposition |4.2.2| Thus, all inequalities are 
equalities, and xs is an optimal pattern. Furthermore, by Proposition 4.2.2 for 


any optimal x there is an optimal E satisfying (3.4|. Note that there are also 


optimal E that do not satisfy (3.41 for any optimal irrigation pattern Xi e. g. it 
is easy to see that E = {x € K" : m^(x) > is such an example. However, 
any optimal E satisfies 

{x S K" : m^{x) > C S C {x G K" : m^(x) > 
for some optimal Xj since for an optimal E the left-hand side and right-hand 


side in (4.11 must coincide and thus all inequalities must be equalities. 


We end this section proving Proposition |3.4.1| 


Proof of Proposition \3.4-l\ The construction of the optimal x from an optimal 
E in the proof of Proposition 4.2.1 is loop-free. Furthermore, it may be chosen 
to have the single path property, since with E fixed, if there are two paths 9i 
and 02 passing through x and y, then without changing the energy O 2 may be 
redirected to have the same path as 61 between x and y. □ 
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4.3 Equivalence of flux- and pattern-based formulation 

Propositions^.3. f] and |4.3.6| together prove 


niin£^’“’^+’^“ [x] = min£’^’“’'^+’'^” [F] 

X T 


as well as the relation (3.3) for two minimisers. 


Proposition 4.3.1. There exists a mass flux T with 


min£:=.“.M+.A«- [^] > ,M- [jT] _ 

X 

Proof. For any ft, > 0 we consider a discrete grid = ftZ" and define discrete 
approximations of via 

ft± = X! + [0: ^]") ^hi , 

iGZ"- 

where 6hi is a Dirac mass centred at hi. Due to the bounded support of ij.±, /i^ 
is a finite weighted sum of Dirac masses. Furthermore, ^ /i± as ft —>■ 0. 

Let X be an optimal irrigation pattern for x^ an optimal 

h h 

irrigation pattern for ^ Further below we will show 

£.£,a,M+.M- > linisupS"’^’'^^^- [x'‘]. 

/i—j-0 

Furthermore, we will later also show that the x^ can be identified with finite 
graphs (or the corresponding fluxes) such that 

Now denoting by Wi the 1-Wasserstein distance, by Remark [3.2.3| we have 
[x^] < aW ijfl) < a/i(j:(K")(2ft + diam(spt/r_|_ Uspt^_)) so that 
the finite graphs have uniformly bounded energy. The corresponding fluxes 
are also uniformly bounded with respect to the total variation norm due to 
||-^G'*llrca < £’p“(G^) and thus are precompact with respect to weak-* conver¬ 
gence. Hence, there is a mass flux T such that ipQh ^ up to a subsequence. 
The lower semicontinuity of the cost then implies 

< liminf ft;p’“(G'*), 

which concludes the proof. 

In order to show > limsup^_,,o [x^], we associate with 

every point x G K" the Lipschitz path I —)■ M" given by 

0 ^{t) = a; -|- t(ft[a;/ftj — x ), 

which connects x with its corresponding point in (here, [cj = ([cij,..., [c„J )^, 
where [c^J = maxjz G Z : z < c^} is the integer part). Now we can define new 
irrigation patterns according to 

fx+(P,3i) iftG[0, |], 
X+(P,-) = X-(g-) = ^'xp(i)’ x'^ipfl) = lx{P,^t-l) iftG(i,|], 

[x-(p,3t-2) iftG(|,l], 
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where i,(t) = 1 — t. It can easily be checked that transports to with 
cost 

= £;;’“(x+) + 5^’“’'^+’'^- [x] + E^p^^ix-). 

The estimate Ep°‘{x±) < ^±{M.'^)ahy/n (since all paths in x± have length no 
longer than hy/n) as well as yield the de¬ 

sired result. 

Finally, we need to show that the x^ can be identified with finite graphs. 
For i,j e Z” let F^ = {p g F : Xp(0) = hi, Xp(l) = hj}. We have (potentially 
changing x^ on a Pr-null set, which does not alter its cost) 

r= U r*,’ 

i.tez’* 


where only finitely many, say N, terms of this union are nonempty, since /i([L 
and consist of only finitely many weighted Dirac measure s. Sin ce x^ 
be assumed to have the single path property (see Proposition3.4.11, x^ ^ 
Lip(/) may be taken constant on each nonempty Fjj, i. e. X^i^ij) = Xij for some 
Xij G Lip(.f)- Furthermore, due to the single path property, the intersection of 
any two fibres Xij{I) ^cid Xki{I) must be connected and can be assigned an 
orientation according to the irrigation direction. Now define for any subset 
S' C Z” X Z" the fibre intersection 


fs= C] U 

(*j)6S ii,MS 


where for simplicity we set Xij{I) = ^ for Fy = 0. There are at most 2^ 
nonempty such intersections fs, and each of them can have at most N connec¬ 
ted components fg,..., fg (again setting some of the fg to the empty set if 
necessary). We have 

x(r,/)= U U /' 

ScZ"xZ" 0<1<N 

with at most N2^ terms being nonempty. Each of the fg can be assigned 
an orientation and a weight Wg = Pp (U(i j)cS f"b) > amount of particles 

travelling on fg (which is constant all along fg). Furthermore, fg must be a 
straight line segment, since otherwise, by straightening the fibres the cost of 
the irrigation pattern is reduced. Hence, we can define a finite graph whose 
oriented edges are the fg, whose vertices are the edge end points, and whose edge 
weights are the Wg. It is now straightforward to check = Ep°‘{G^) 

as required. □ 

The proof of the opposite inequality requires a few preparatory lemmas. 

Lemma 4.3.2 (Almost a F-convergence lemma). Suppose that 

• /^+, P-- are discrete measures such that ^ p- as N ^ oo; 

N N 

• En is a minimiser of ; 
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Then, !F is a minimiser of . 

Proof. To achieve a contradiction suppose that J- is not a minimiser of , 

that is there exist iF' such that div = /r+—and (J"') < a'+.m- 

On the right-hand side, we have (up to a subsequence) 

< liminf£:^’“’^?’^~(J'Ar) 

N^oo 

due to the weak-* lower semicontinuity of the energy. Thus, given ry > 0, for 
large N we have 

On the left-hand side, by definition of ifp^, there exists a sequence such 
that 

N N 

lim (J-') -h ? 7 . 

Af—>00 

Choosing ry suitably close to 0 now yields {iFj^) for 

N large enough, a contradiction to the optimality of Tn- □ 

Definition 4.3.3 (Paths in a graph and their weight). Let G be an oriented 
graph with edge set E{G). A path in G is a sequence ^ = (ei,..., Cfc) of edges 
ei,..., Cfc G E{G) such that e~ = ef_i for i = 2,... ,k, where e~ and e"*" denote 
the initial and final point of edge e. 

Suppose that G is also weighted and has no cycles. Let us denote by S(G) 
the set of maximal paths on G, that is, paths that are not a subsequence of any 
other path. 

The weights w{£,) of all paths ^ G S(G) is defined by the system of equations 
li'(e) = w(^), e G E(G), 

eG? 

whose solvability follows from |Xia031 Lemma 7.1]. 

Finally, for Sq C S(G) we define 

|So| = w(f). 

«6Ho 


Lemma 4.3.4 (Bound on fibre length). Let /i+,/r_ G fbm(K") be discrete 
measures of equal mass with support in a convex set £l C K" and let G be any 
discrete mass flux between and yt_. Let Sq denote the set of paths in G of 
length greater than 2adiamn, and let G(So) denote the graph whose associated 
vectorial measure is given by 


-^G(Ho) = ^w(e)Me 
563o e6j 


(cf. Definition 2.1.,^ . Then there exists a discrete mass flux G' between /r_|_ and 
pi_ such that all its paths have length bounded by 2adiamf2 and 


Gp’“(G)-Ap’“(G') > ||J^G(Ho)l|rca “ adiamn|.=o| > adiamfl|So|, (4.2) 

II^G - -AG'lIrca < II-Ag(Ho) lirca + adiamfl|So| < 3(G;’“(G)-£;;’“(G')) . (4.3) 
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Proof. Definitions |2.1.3| and |2.1.4| can be extended to paths and thus allow to 
define graphs as sums of edges and of paths, which we make use of in the 
following to simplify the exposition. Let 

G'= Y. MOi + G"= Y «^'(e)e + G", 

CGH(G)\Ho eGE(G) 


where 


w'{e) = w{e) - Y 

e6«6Ho 


and G" is a graph composed of straight edges, which recovers the flux conser¬ 
vation condition div G' — div G = /r+ — /r_. Denoting the length of an edge e 
by l{e), we can now compute 


£;|’“(G) - E\ 


(GO 


eeE(G) 


= E 

eGE{G) 


Y 

eG?GHo 


lie) - l;^’“(g") 


«eHo eeC 


> ll-^GCHoillrca “ adiamSl|So|. 


The relation 


lie) > 2adiamf2|So| now con- 


^eG£ 


eludes the proof of (4.21. Equation (4.3 \ directly follows from G—G' = Gi^^o) — 
G". □ 


Finally we will need the following compactness lemma for transport path 
measures. 

Lemma 4.3.5 (Compactness for transport path measures). Let G > 0 and 
LI C M" be compact, and consider the set 

©c = {0 e 0 : 0(7) C n and |0(f)|dt < g} C 0. 

Let rjN G TPM(/i+,/i_) be a sequence of transport path measures such that 

r?7v(©\ ©c) = 0. 

Then, up to a subsequence, rjpf ^ rj in the sense 


(fiO) d??Ar(0) 


(pi0)drii0) 


for all (fi G Cbi &), 


where C(,(©) denotes the set of bounded continuous functions on ©. Moreover, 
V G TPM(/r+,/x_). 

Proof. Note that © is separable (which follows from the separability of Lip(7)) 
and that ©c is a (sequentially) compact subset of ©. Indeed, let 0„, n = 1, 2,..., 
be a sequence in ©c. Upon reparameterisation of each element (which does not 
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change the sequence), the are uniformly Lipschitz. Thus, by the Ascoli- 
Arzela Theorem, up to a subsequence we have t 0 S C(J; 17). Furthermore, 


f | 0 (t)| dt < liminf f | 0 Ar(t)| dt < C. 

Ji ^ JI 

As a consequence, the are all supported on the same compact set and are 
thus tight (i.e. for every e > 0 there exists a compact such that ri]\[{K^) < e). 
Furthermore, due to rjN G TPM(/j,+ ,/j,_) they all have the same mass. Hence, 
by Prokhorov’s Theorem (which assures weak compactness for a tight set of 
measures; see |Bil99| l we get ^ ?7 up to a subsequence, as desired. 

It remains to prove po^rj = /i_|_ (the proof oipi^r] = works analogously). 
Because of Po^rjN = P+ for all N we have 


[ <p(Po(^)) d? 7 Ar( 6 ») = [ (p{x) dp.+ {x) for all ip G Cb{n) . 

Je Jn 

Due to ryjv ^ ?7 as well as tpo pQ G Ct,(0), letting TV —>■ oo we finally arrive at 

[ ip{po{9)) dr]{9) = [ ip{x)dfi+{x) for all ip G C(,(n), 

Je Jn 


that is, po^p = P+- n 

Proposition 4.3.6. We have 


min^e.a.M+.M- r^i < niin£:^’“’^+’'"- m . (4.4) 

X 

Furthermore, for any optimal mass flux T there is an optimal irrigation pattern 
X so that both are related via 




dF = 


J J <p{Xp{t)) ■ Xp{t)dtdPr{p) for allip (4.5) 


Proof. In the first part of the proof, we construct a pattern x from an optimal 
flux. So let the flux T be optimal. We may assume ^ since 

otherwise there is nothing to show by Remark |3.2.3| Let a sequence of finite 
weighted graphs such that 

• Gtv is a discrete mass flux between some and 

• {p^,p^,Pgn) ^ 

• ^^’“(GAr) 

^ N N 

Note that if Gat or rather Pq^ is a minimiser of , by Lemma 

must have 

lim A^’“(Gjv) - EI.'^{Gn) = 0 . 

N—^oo 

Let Gn denote the cycle-reduced graph Gat. By Lemma|4.1.9| ifp“(GAr) < 
F;^’“(GAr) and 


4.3.2 


we 


II-^Gn 


Pa 


< E^f\Gn) - E^^Gn) < E^^%Gn) 


||rca _ -^p 


E^'^Gm) 


0 . 

N—¥00 
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Thus, without loss of generality, we may replace the Gn by discrete mass fluxes 
without cycles, and from now on Gjv is supposed to be without cycles. 

We may even assume the Gw to only contain paths with length bounded by 
2adiam(n), where is a ball containing spt^+ and spt/i- . Indeed, let G^ and 
be the graph and the set of paths from Lemma 4.3.4 associated with Gw. 
By Lemma [4. 3. 4| we obtain 




< 3(G;’“(Gw)-i;r(GW)) < ?.{E^^\Gn)-E^^\Gn)) 


N—^oo 


0 . 


Thus we may replace the Gw by the G'^, which have uniformly bounded path 
lengths. 

Summarising, from now on we may assume the Gw to have no cycles and to 
have path lengths bounded by 2adiamn. Hence, by Remark [4.1.10| there exist 
corresponding transport path measures ?7w G Since they para¬ 

metrise the graphs Gw, the 77 W have support on paths with lengths bounded by 


2adiamS2 and images in i32adiama(^)- Thanks to Lemma 4.3.5 we thus have (up 
to a subsequence) ijn ^ rj for some t] G TPM(/j,+ , p,_). By Skorokhod’s Conver¬ 
gence Theorem |BCM09l App. A, Thm. A. 8 ], there exist a sequence of irrigation 
patterns xn parameterising r]ff and an irrigation pattern x parameterising 77 
such that 

Xn{pg)^-^ x{pg) for all per. 

By Proposition |3.2.8] we have 

< liminf £;;’“(xw) = liminf £;^’“(Gw) = [E ], 

N —¥00 N —¥00 


and (4.41 is established. 


In the second part of the proof we now explain the relation given by formula 


(4.5) between the constructed x E. After reparameterisation according to 


Proposition3.2.5 we may assume the xn{p,‘) to be uniformly Lipschitz with 
constant 2adiam(r2). Thus, for each p we may extract a subsequence converging 
weakly-* in against x(p, •)• Since any subsequence contains such a 

converging subsequence, actually the whole sequence xn{p, ■) converges weakly- 
* against x{Pi ’)• Now for any (p G Cc(K";M") we have (the second equality can 
be easily verified edge by edge) 

/ (p-dJ'= lim / p • dJ 'gjv = Jim / / (p(xn(p, t)) ■ Xn(p, i) dt dPr(p) ■ 

Note that p{xn{p,‘)) converges in while xn{p,‘) converges weakly-* in 

L°°{I) so that 

p{xn{p, t)) ■ Xn{p, t) dt =: Jw(p) J (p) ■= J p{xiP, t)) ' x{P, t) dt . 

Together with the uniform bound Jw(p) < \\'p\\l°°'E^{xn{p, I)) < 2||(p|jiooadiam(Sl), 
this allows application of Lebesgue’s dominated convergence theorem, from 
which we Anally obtain 

J^(p-dE = \im^ J J^(p{xN{p,t))-XN{p,t)dtdPr{p) = J J p{x{pE))-x{pE) dtdPr{p), 


the desired formula. 


□ 
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